Evidence of anisotropic magnetic polarons in Lao.94Sro.o6Mn03 by neutron scattering 

and comparison with Ca-doped manganites 
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Elastic and inelastic neutron scattering experiments have been performed in a Lao.94Sro.o6Mn03 un- 
twinned crystal, which exhibits an antiferromagnetic canted magnetic structure with ferromagnetic 
layers. The elastic small q scattering exhibits a modulation with an anisotropic q-dependence. It 
can be pictured by ferromagnetic inhomogeneities or large magnetic polarons with a platelike shape, 
the largest size («17A) and largest inter-polaron distance (~ 38A) being within the ferromagnetic 
layers. Comparison with observations performed on Ca-doped samples, which show the growth of 
the magnetic polarons with doping, suggests that this growth is faster for the Sr than for the Ca 
substitution. Below the gap of the spin wave branch typical of the AF layered magnetic structure, 
an additional spin wave branch reveals a ferromagnetic and isotropic coupling, already found in 
Ca-doped samples. Its q-dependent intensity, very anisotropic, closely reflects the ferromagnetic 
correlations found for the static clusters. All these results agree with a two-phase electronic seg- 
regation occurring on a very small scale, although some characteristics of a canted state are also 
observed suggesting a weakly inhomogeneous state. 



PACS numbers; 74.50. C, 75.30. K, 25.40.F, 61.12 
I. INTRODUCTION 

In the simple ionic picture of Lai_2;Ba;Mn03 (B=Sr^+, 
Ca^"*", or Ba^+), upon increasing x, the Mn'^^ ions char- 
acterized by three t2o and one e^ electrons are gradu- 
ally replaced by Mn + with empty e^ orbitals. These 
compounds have attracted much attention due to their 
anomalous transport, magnetic and structural proper- 
ties. The low-doping range of these compounds preceding 
the ferromagnetic and metallic phase is far from being 
fully explored, and the mechanism driving the metal- 
insulator transition is still unclear. On the theoretical 
side, the situation is rather controversial. Two pijCtures 
have been proposed. In the model of de Genncad, the 
competition between the AF superexchange and the fer- 
romagnetic double-exchange couplings associated with 
the hole hoppingj, results in a canted state. In con- 
trast with this picture, several authors have found that an 
inhomogeneous representation of ferromagnetic hole-rich 
regions (called ferrous or ferromagnetic polarons) , within 
an antiferromagnetic hole-poor mediiim, could be more 
favorable than that of a canted statdjQ. Similar reaulte 
have been obtained by Monte-Carlo simulationsfliliJ. 
This picture of ferromagnetic polaron, here, differs from 
the spin polaron picture discussed by Kasuyafl, which 
consists of a ferromagnetic cloud surrounding one hole. 

In Lai_xCaa;Mn03, interesting observations have been 
recently reported: a "static" ferromagnetic modulation 
has been observed at small wave- vectors q, interpreted 
in terms of magnetic inhomogeneities or droplets, with 
a distribution analogous to that of a liquidM. More- 
over, below the gap of the spin wave branch with an 
anisotropic dispersion typical of the AF layered struc- 
ture, an additional spin_sY*5^ branch with an isotropic 
dispersion was observedlij113. Whether these findings 



are restricted to Ca substitution, in relation with the 
stronger induced lattice distortion, or are more general, 
is an important question. The present study carried out 
on Sr-substituted samples sheds light on this issue. Sr 
is known to induce a smaller distortion of the Mn-O-Mn 
bonds than Ca, resulting in a higher hopping probabil- 
ity, t, for the Cg holes. Therefore, in the ferromagnetic 
and metallic state, at large x values, the ferromagnetic 
transition temperature is much higher for Sr than for 
Ca substitution. This enhanced ferromagnetic character 
can also be seen at smaller x, by a shift of the critical 
concentration Xc for the metal-insulator tcansition. For 
instance, the X5r=0.175 sample is rrvatallicoEj, whereas 
the xca=0.194 one is still insulatingla. Previous inelas- 
tic neutron iiaeasurements-pn low doped Lai_a:Sra;Mn03, 
with x=:0.05llZl and 0.09t2l, have reported a single spin 
wave branch, with AF interplane and F in-plane mag- 
netic couplings. 

The present study, carried out on an untwinned 
xsr=0.06 sample, reveals two additional features. Below 
the temperature of the magnetic transition, a ferromag- 
netic modulation of the intensity of neutron scattering 
can be observed at small q, similarly to pur previous 
observations in one Ca-substituted sampleEJ. The un- 
twinned sample reveals the anisotropic character of this 
modulation. The size of the inhomogeneities is about 
twice as larger (17^1) within the ferromagnetic layers as 
perpendicularly to them, indicating platelets instead of 
droplets. Furthermore, an additional spin wave branch 
is observed with an isotropic dispersion curve, charac- 
teristic of a ferromagnetic coupling. The q dependence 
of the corresponding susceptibility, along the two main 
symmetry directions, reflects closely the correlations of 
the static inhomogeneities, showing the same anisotropy. 
This study confirms the general character of our previous 



observations. The density of the platelets which is much 
lower than the hole density confirms a picture of an elec- 
tronic phase segregation on the scale of two tens of A. 
Comparison with the previously studied Ca-substituted 
sample (xpa^O.OS) and with a new one with a smaller 
xca value (xca=0.05) reported here , brings up to de- 
duce that the fcrromagnctism induced by hole doping is 
enhanced upon Sr substitution. Such an enhancement is 
observed both in the size of the ferromagnetic polarons 
and in the characteristics of the spin dynamics. 

The rest of the paper is organized as follows: In section 
II the magnetic structure is determined. In section III A 
and B, elastic diffuse scattering data at small angle and 
around the (110) Bragg peak are reported. In section IV, 
the spin dynamics is investigated. Finally, in section V, 
the overall results are discussed. 



II. MAGNETIC STRUCTURE 

Experiments have been carried out on the three-axis 
spectrometer 4F1, installed on a cold source at the re- 
actor Orphee (Laboratoire Leon Brillouin), on an un- 
twinned single crystal of Lao.94Sro.o6Mn03 (0.1 cm'^), 
grown by a floatting zone method. The nominal x=0.06 
value agrees witlt-^ deduced from the lattice param- 
eters reported int^EJ. The structure is orthorhombic 
(c/v2 < a < b) with Pbnm symmetry. The temperature 
dependence of (110), (001) and (112) Bragg-peak inten- 
sities was studied. The r=(001) Bragg-peak intensity 
determines a long range antiferromagnetic order below 
Tat^ISOK (see Fig. 1-a). The evolution of the r= (112) 
Bragg peak intensity, very sensitive to the occurence of 
ferromagnetism because of its small nuclear intensity, de- 
termines a long range ferromagnetic order at T/f =124 K, 
smaller than T^r (see Fig 1-b, c). The absence of any 
magnetic component below Tk in the (002) Bragg peak 
intensity, indicates that the ferromagnetic spin compo- 
nent is along [001] (geometrical factor). Therefore an 
average canted state is defined below Tk , which consists 
of ferromagnetic layers in the (a, b) plane, antiferromag- 
netically stacked along c, and with a weak ferromagnetic 
spin component along c as sketched in Fig. 1-d. The 
two transition temperatures agree with the phase dia- 
gram of Ca-doped samplesEj being between those found 
for xca=0.05 and xca=0.08. However, an interesting dif- 
ference between the two systems appears for the spin 
deviation or canting angle 9 deduced from the weak fer- 
romagnetic component. From the (112) Bragg peak in- 
tensity, an average deviation angle 9 from b towards c 
of 13° is determined at T=15K; this value is the same as 
that found for xca=0.08, and much larger than the value, 
9—5'^, determined for xca=0.05. Therefore, the increase 
of the canted angle 9, upon increasing x, is faster for a 
Sr substitution than for a Ca one. It is worth-noting 
that the Bragg-peak intensities define an average mag- 
netization TO but do not inform about its homogeneous 



or inhomogeneous character. The inhomogeneities con- 
tribute to diffuse scattering between Bragg peaks which 
is studied below. 
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FIG. 1. a) Temperature variation of the intensity of (a) 
the (001) Bragg peak defining Tiv (130K) , and (b, c) the 
(112) Bragg peak defining Tk= (124K). In (d), the drawing 
represents the average magnetic structure. 



III. ELASTIC DIFFUSE SCATTERING 

Elastic (?iti>=0) diffuse scattering measurements, which 
give access to short-range ordered correlations, have been 
performed at several temperatures along several Q-paths 
in the reciprocal lattice, with Q=q-|-T. In section A, dif- 
fuse scattering intensity close to the direct beam t=(000) 
is reported along the main symmetry directions. In sec- 
tion B, diffuse scattering around x=(110) is reported 
along the same q directions. 

These two types of experiments appear complemen- 
tary. The scattering intensity close to the direct beam 
is more easily modelized, but, because of the geometrical 
factor, the [001] symmetry direction is unobservable. In 
contrast, this direction can be observed close to t=(110). 
We also compare the diffuse scattering close to t=(110) 
with that measured close to t=(220), where, due to the 
Mn magnetic form factor, the magnetic intensity is ex- 
pected to decrease. Finally, antiferromagnetic correla- 
tions have also been measured close to t=(001) along 
the [001] direction and reported in section C. 



Ferromagnetic short-range order studied close to 
r=(000) 
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FIG. 2. Elastic magnetic scattering Ii5x-Iht versus q 
in barns per Mn (left hand scale) and neutron counts (right 
hand scale): (a) along [110] and [112] q directions and (b) 
along [001] q direction. The lines are guides to the eye. (c): 
schematic representation of a modulation of longitudinal com- 
ponents S^. 

Small angle scattering experiments have been per- 
formed in the range 0.lA~^<q<0.6j4~^ , and in the tem- 
perature range 15K<T<300K along [110], [112] and [GDI] 
directions, using ki—1.25A~^. 

As inEJ, the magnetic intensity I(q) was obtained 
by subtracting the high temperature spectrum, 1ht{<i), 
found to be q-isotropic and temperature independent for 
T>Tx, from the spectrum measured at 15K. Such a 
subtraction procedure assumes that lHT{<i) has a nu- 
clear origin assigned to dislocations. Actually a resid- 
ual scattering intensity above Tk can be observed close 
to x=(110) where dislocations do not contribute (sub- 
section B). Since it is very small, the subtraction proce- 



dure is valid. I(q) was converted in barns per Mn us- 
ing a Vanadium standard, and is reported in Fig 2-a for 
q along [110] and [112] and in Fig 2-b along [001]. In 
both [110] and [112] directions, a modulation is observed 
aromid q=0.2A~^, differing in intensity. The weak maxi- 
mum observed at larger q along [110], since not observed 
close to the equivalent symmetry point (110), is likely 
a spurious effect. The magnetic contribution observed 
along [001], is nearly zero around 0.2A~^, and increases 
shghtly with q (Fig 1-b). 

2) Analysis r—. 

In our previous study of a Ca-doped sampldl^, the ap- 
parent isotropy was the consequence of twinned domains. 
We have proposed a crude model of magnetic inhomo- 
geneities or droplets. The basic assumption is that the 
scattering intensity results from the contribution of lon- 
gitudinal spin components S^ along [001] as sketched in 
Fig. 2-c. In addition, we have assumed that the char- 
acteristics of the magnetic inhomogeneities (shape, dis- 
tance) are the same whatever the direction (isotropic 
model). Our untwinned sample provides the opportu- 
nity to check these two assumptions. By analogy with 
the chemical problem where the nuclear scattering length 
hi J replaces Sf.,-, the Fourier transform I(Q) of the cor- 
relations of S^ spin components may be written: 



E::,=i(5^5f 



(S"'-^)2)e*'5(^-"^j) 
Here, the symbol _L means the projection of S^ in the 
plane perpendicular to Q. The first term gives rise to 
the ferromagnetic Bragg peaks and the second one to 
ferromagnetic short range order. The bar stands for spa- 
tial average. The thermal average has been omitted. In 
the present case, Q=q, since T=(000). By approximat- 
ing the inhomogeneous character of the magnetization to 
a step function characterized by two magnetization val- 
ues (see Fig 2-c), and replacing the discrete spin Si by 
the continuous magnetization density m;^(r), the second 
term may be factorized into two functions F{q, R) and 



J(g, d, 



m-) ^min 
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I{q) = CT^{Q)\F{q,R)\^J{q,d„„d„,,r.) 



(1) 



where F{q, R) is the Fourier transform of one clus- 
ter assimilated to a sphere of radiua_j?, and J(g) de- 
scribes the cluster spatial distributiontiJ. Its full expres- 
sion which uses the average inter-cluster distance dm and 
a minimal jdistance of approach dmim has been previ- 
ously givenlld. ^{Q) is the Mn magnetic form factor. 
Close to T=(000), T{Q) « 1. From equation (1), mainly 
two length scales are derived: (i) the typical size of the 
inhomogeneity R is mainly deduced from the q depen- 
dence beyond the maximum of I(q), since J{q) in (1) 
tends asymptotically to unity, (ii) the characteristic dis- 
tance between the inhomogeneities, dm, (or equivalently, 
the inverse of the density Ny^), deduced from the liq- 
uidlike distribution function J{q,dm,dmin) in (!)■ An 
approximate determination can also be made from qmaa; 
(q at maximum of intensity) with dm='^T^ l'\max- These 



two parameters are included in the scale factor C of (1), 
C = (7ro)^|ATO^"'"pAfyV^, (7, ro being constants for neu- 
tron scattering). The only free parameter when fitting 
is the squared magnetic contrast |Am^-'"p, given by the 
difference of the magnetizations between the two regions 
(sketched in Fig. 2-c). In this model, the intensity is 
very sensitive to the size (R enters at the power 6 in the 
constant C for a spherical shape through the squared vol- 
ume V'^). Therefore, in an isotropic model, the scatter- 
ing intensity I(q) changes only by the geometrical factor 
as the q direction varies, decreasing by «l/2 from [110] 
(perpendicular to c) to [112], down to zero along [001]. 

Comparison of these predictions with observations of 
Fig. 2-a and Fig. 2-b, leads to the following conclusions. 
The striking result is the absence of any intensity in the 
range 0.1<q<0.2 A~^ along [001], where the maximum 
of intensity was observed in twinned samples. This con- 
firms the assignment of this magnetic scattering to S^ 
spin components. The increasing intensity at larger q 
indicates short-range antiferromagnetic correlations (the 
AF Bragg peak (001) is located at q=0.83A-^ in Fig 2). 
They will be discussed in section C. Along [110] and [112], 
the scattering intensities do not only differ by a simple 
geometrical factor, but also by some anisotropy in their q 
dependence. The slower q dependence of I(q) beyond the 
maximum, for q along [112] indicates a smaller size along 
the [112] direction (at « 45° from c) suggesting a still 
smaller size along the [001] direction. Diffuse scattering 
experiments close to t=(110) reported in section B will 
characterize this anisotropy. 

In the absence of any anisotropic model, a semi- 
quantitative description can be made, using the equa- 
tion (1) to fit I(q) along the [110] direction only. The fit 
(continuous line in Fig 2-a) provides a correlation length 
2R—^//fiil7A and a density of clusters Ny which cor- 
responds to an inter-cluster distance which is twice the 
diameter value. The low density of the ferromagnetic 
clusters cannot be explained by the picture proposed for 
diluted chalcogenide alloyaJ where a single electron fer- 
romagnetically polarizes the neighbor spins but rather 
indicates an electronic phase separation. In the rest of 
the paper, these ferromagnetic inhomogeneities are called 
ferromagnetic clusters or large ferromagnetic polarons in 
an equivalent way. 

In Fig. |-3. I(q) for X5,.=0.06 is compared with I(q) for 
xca=0.08liJ and xca=0.05, in absolute units, for q along 
the [110] direction and at T=15K. For the Ca-substituted 
samples, a factor 2/3 takes into account the fact that 
one domain out of three does not contcibute (e.g. we 
report 3/2 I(q), I(q) being displayed inO). We estimate 
a possible error of 10% in the intensities, due to the fact 
that they are corrected for sample transmission which is 
calculated instead of measured. 

Considering first the two Ca-substituted samples, the 
analysis (dashed and dot-dashed lines) shows the growth 
of the clusters with xca from « 14 A to « 17 A as xca 
varies from 0.05 to 0.08. In addition, there is a small 
shift of <imax with X(7a, interpreted as an increase of 



the cluster density. The increase of the intensity by a 
factor two nearly agrees with equation (1) , being 30% 
smaller than expected. Interestingly, when comparing 
the xsr—O-QG sample to the two Ca-substituted ones, the 
maximum intensity of I(q) (« 2 barns/Mn) and the q de- 
pendence beyond qmax nearly coincide with the values of 
the X(7a=0.08 case. Within the above model (continuous 
and dot-dashed lines), it indicates that the cluster size is 
the same in the two samples. We note a small shift of the 
Imax value between both samples, also observed in the 
raw data around t=(110), reported in the inset of Fig. 
3 in reduced lattice units (rlu), which could indicate a 
slightly smaller cluster density in the xg^^O.OG sample. 

It results from the analysis of these three samples, 
that the intensity seems to be mainly determined by 
the size, in agreement with the model given by equa- 
tion (1). This also indicates that the magnetic contrast 
value is Ihe same in the three samples. As previously 
reportcdiiij, the value deduced from the scale factor C 
in (1), |Am^^|R:!0.7/Li_B is surprisingly small. This could 
result from an overestimate of the volume due to the 
shape anisotropy, the [110] direction corresponding to the 
largest size. Moreover, a slow decrease of the magneti- 
zation starting from the core of the polaron, as sketched 
in Fig. 2-c, is more realistic than the step function of 
the present model. We indicate that NMR experiments 
under an applied field suggest that some spins are fully 
ferromagneticLJ. 

In conclusion of the above analysis, the size of the clus- 
ter, at a given x, depends on the chemical dopant, being 
larger for Sr than for Ca. This outlines the role of the 
hopping probability "t", regulated by the mean Mn-0- 
Mn bond angle, to define the size. 
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FIG. 3. Magnetic scattering intensity versus q in 
barns/Mn. Full circles: xsr=0.06, open circles: xca~0.08 
and open squares: xca=0.05. The continuous, dot-dashed 
and dashed lines are fits with the model (see the text). Inset: 
diffuse scattering (raw data) close to the (110) Bragg peak in 
equivalent reduced lattice units (riu). Full circles: X5r=0.06, 
open circles: xca=0.08. The continuous lines are guides for 
the eyes. 



B. Ferromagnetic short-range order studied close to 
r=(110) 

Diffuse scattering intensity has been measured close to 
T=(110) for main symmetry directions and as a function 
of temperature, using the incident neutron wavevector 
k,=l. 551^1. 



FIG. 4. a: Raw intensity versus q, converted in A . q in 
rlu has been deduced from q=Q-T, t=(110). a: along [110], 
(0.4,0.4,0)<Q<(1,1,0). b: along [001], (1,1,0)<Q<(1, 1,0.85) 
andc: along [112], (1,1,0)<Q< (1.3,1.3,0.6). The continuous 
lines are fits with a gaussian (modulation ) and delta func- 
tions (Bragg peak) except at 130K in (a) where a lorentzian 
function is used. In b, the gaussian component, corrected for 
the geometrical factor, is shown by a dotted curve. 
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In Fig. 4 a-c, raw data are reported as a function of 
q in A-i, for q//[110] (a), q//[001] (b) and q//[112] (c) 
at several temperatures. Here q is defined with respect 
to the (110) Bragg peak. Along [110], a modulation of 
I(q) appears very well-defined at 16K. It agrees with the 
observations reported close to t=(000) (Fig. 2), consid- 
ering that slight differences may be expected due to the 
subtraction procedure. Along [112] and [001] a modu- 
lation of I(q) is also clearly seen at 16K, growing below 
the transition temperature, but with a smaller amplitude. 
Note the change of scale in both coordinates between Fig. 
4-a and Fig. 4-b. Along [001], beyond q=0.55A"\ the 
increasing intensity corresponds to the tail of the short- 
range antiferromagnetic order as q approaches the AF 
(111) Bragg-peak. 

By neglecting the displacement of Mn atoms, we can 
use the same analysis as that performed in section A. 
Since S*^ in equation (1) is the spin component per- 
pendicular to Q— q+T with t=(110), all directions are 
observable. Moreover, the corrections due to the geomet- 
rical factor can be neglected (cf Fig. 4-b). 

All features indicate that the ferromagnetic corre- 
lations which define the ferromagnetic cluster have a 
smaller extension along [001] (antiferromagnetic stacking 
of the layers) than along [110] (within the ferromagnetic 
layer). Assuming that all q directions within the ferro- 
magnetic layers are equivalent, such an anisotropy defines 
a plate-like shape for the inhomogeneities or that of an 
ellipsoid of revolution with main axes ^//, S,// , ^_l- Since 
few directions have been studied, the q-dependences have 
been fitted independently with a gaussian function as an 
approximation of the F(q,R) function used in section A. 
Along [001], the q dependence beyond qmax determines a 
typical size £^±k7A (« 2 Mn-Mn distances), nearly twice 
as small as the value ^//«17A found along [110]. The 
origin of this anisotropy will be discussed below, in rela- 
tion with the anisotropic double exchange due to the AF 
layered structure. 

In addition, the qmax value at the maximum of inten- 
sity is larger along [001] (q« 0.25^1^^) than along [110] 
(ss G.ISA^'^). We deduce an approximate inter-cluster 
distance dm—^TT / q„iax smaller along [001] (« 25A), than 
along [110] («38l). 

The behaviour with temperature has been mainly stud- 
ied for q along [110] (see Fig. 4-a). As T increases, a 
decrease of the intensity is observed. 

The decrease of the elastic scattering does not result 
from a dynamical broadening. This was checked by en- 
ergy scans at one constant q value, using the energy 



resolution 40/xeV (using an incident neutron wavevector 
ki=l.l A~^). Therefore, the signal remains "elastic" , 
even for our best energy resolution up to Tk- Close to 
Tk a large scattering intensity enters in our finite experi- 
mental window, as shown in Fig 4-a at T=130K. It can be 
fitted by a lorentzian function, and is attributed to crit- 
ical scattering. Above Tk a residual bump of intensity 
persists, centered at q=0. In conclusion, the scattering 
intensity attributed to ferromagnetic polarons vanishes 
by approaching Tk as does the ferromagnetic Bragg-peak 
intensity. The evolution of the q-dependence with tem- 
perature is difficult to determine, except at small q where 
a shift of qmax to q=0 is observed. The residual scatter- 
ing above Tk suggests that some electronic segregation, 
which, in the temperature range T<Tk gives rise to the 
ferromagnetic clusters, could also persist above Tk- Ac- 
tually, such a persistence can be better studied from the 
temperature evolution of the unusual spin choiamics, as 
was reported for the Ca substituted samplealj. 
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FIG. 5. Comparison between diffuse scattering observed 
along [110], around t=(110) (empty circles) and (220) (full 
circles) using fci=2.662A~^. Close to r=(110), only a modu- 
lation fitted by a Gaussian (continuous line) is observed with 
the Bragg peak (dotted line). Close to r=(220), a new com- 
ponent (dotted dashed line) is observed in addition to the 
magnetic modulation (dashed line). 

Finally, in Fig 5, the diffuse scattering along [110] 
observed around t=(110) and r=(220) are compared, 
both being measured in the same experimental condi- 
tions {ki=2.662A~^). The modulation does not readily 
appear around t=(220). Actually, besides the contribu- 
tion of the Bragg peak, I(q) consists of a q-modulated in- 
tensity, twice smaller than that observed around r=:(110) 
(dashed line) , superimposed to a huge intensity at smaller 
q, approximated by a lorentzian function (dotted -dashed 
line). The decrease of the modulated intensity with Q 
as T increases from (110) to (220), agrees with the de- 
crease of the squared Mn '^+ magnetic form factor, Tifff' 
in equation (1), emphasizing the magnetic origin of the 
modulation. The lorentzian contribution, since not ob- 



served around r=(110), has no magnetic origin. It reveals 
a deviation from the perfect atomic long range order. 
Deviations from atomic positions give rise to a nuclear 
intensity varying path Q^, whatever their isolated or seg- 
regated charactered. We recall that diffraction measure- 
ments performed on a very large Q range, and using a 
pair distribution function have determined two distinct 
Mn-0 distances attributed to Mn''+ and Mn^+N, and 
therefore a disorder for the oxygen positions. The present 
observation could be related to this disorder. 



C. Antiferromagnetic short-range order 

As a feed-back effect of the ferromagnetic inhomo- 
geneities described above, we also expect to observe some 
AF short-range correlations within the canted AF layered 
structure. Elastic diffuse scattering has been measured 
along [001], close to the antiferromagnetic Bragg-peak 
T=(001) at T=16K and T>TAr (T=200K) and reported 
in Fig. 6. There, due to the geometrical factor, only 
correlations of spin components transverse to the [001] 
direction (S^^, S^) are measured. Actually, below Tk, 
a diffuse scattering intensity grows. At 16K, the inten- 
sity consists of a weak modulation in the same q range 
as in Fig. 4-b, superimposed to a q-lorentzian intensity 
at the foot of the Bragg-peak. This latter reveals an 
additional magnetic disorder, corresponding to a weak 
deviation from the AF ordering on a scale larger than 
the ferromagnetic cluster (wl30A from the q-linewidth) . 
It likely corresponds to some rotation of the transverse 
components in the (a, b) plane around their mean b di- 
rection, which is an additional degree of freedom for these 
spin components. 
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FIG. 6. Scattering intensity versus q in A~^ along [001], 
with q=Q-r, t=(001): filled circle: T=16K empty circle: 
T=200K. Data at 16K are fitted by the sum of a delta func- 
tion (dotted curve) for the Bragg peak with a lorentzian 
(long-dashed line) and a gaussian (small-dashed line) func- 
tions. 



IV. SPIN WAVES 

In the previous section, we have described an inho- 
niogeneous state from a phenomenological point of view: 
cluster size, shape anisotropy etc. The determination of 
the magnetic excitations allows a better understanding of 
this peculiar magnetic state. Moreover, in the absence of 
twinned domains, the magnon intensities can be precisely 
determined, bringing new important informations. Since 
we have already studied them in the Ca-doped samples, 
we can get a quantitative comparison between the two 
systems. 
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FIG. 7. Dispersion curves of spin waves for q along [001] 
and along [110]. Filled symbols correspond to the high energy 
spin wave branch, and empty ones to the low energy one. The 
continuous lines are guides for the eye. The arrows point the 
gap value of the low-energy branch for xca=0.05 (upper one) 
and xca=0.08 (lower one). 

Inelastic neutron scattering experiments were per- 
formed on the spectrometer 4F1, at constant final neu- 
tron wavevectors kf — 1.55A^^ and kf = 1.3A^^ at 50K 
and 15K. Energy spectra for excitations in creation mode 
process have been fitted by a lorentzian function as spec- 
tral function, convoluted with the spectrometer resolu- 
tion: 
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From the fits, the energy modes 'jj(q), their damping 
F(q) and the energy integrated intensity x(Q) are deter- 
mined. This latter parameter contains the geometrical 
form factor and any Q dependence related to the sym- 
metry of the magnetic structure. Since the modes are 
underdamped, the choice of the spectral function does 
not affect the results. n{'huj) is the Bose population fac- 
tor. 

To benefit from the increase of intensity related to the 
thermal population factor, the two main-symmetry di- 
rections [001] and [110] have been fully studied at 50K, 
whereas only the [001] direction has been studied at 15K. 



In Fig. 7, huj{q) are reported for q in rlu along [001] 
and [110] directions at T=50K, focusing on the energy 
range below 3 meV. The scale in abscissa has been cho- 
sen to allow a direct comparison of the two directions 
in A~^ units. The similarities of the dispersion curves 
with previous observations in Ca-substituted samples are 
strikingi3il3: two dispersive spin-wave branches are ob- 
served. One with a large gap and an anisotropic disper- 
sion, characteristic of the AF layered structure. Cor- 
responding spectra are shown in Fig 8-a for direction 
[001]. The other one, with a small gap and an isotropic 
dispersion, is observed at nuclear Bragg-peaks only. It is 
characteristic of a ferromagnetic couphng. Along [001], 
it levels up at the zone boundary (001.5), and lies be- 
low the gap of the former (see the hatched area along 
the [001] direction in Fig. 7). Corresponding spectra are 
shown in Fig. 8-b, and 8-c for directions [001] and [110] 
respectively. 

By fitting the high energy branch along [001] with a 
dispersion relation deduced_from a Heisenberg hamilto- 
nian as in the pure systemEj, we determine an interlayer 
"effective" antiferromagnetic coupling (J2=-0.3 meV at 
T=15lO, strongly reduced with respect to the undoped 
systemEa (J2=-0.58 meV at T=15K). Along [110], high 
energy measurements are lacking to determine the Ji 
coupling, which was found to increase slightly with x in 
Ca-doped samplesll3. 

The low-energy spin wave branch has been analyzed 
with a phenomenological quadratic law at small q, 
hLo=Dq^ -^huoQ {D — GJSuq, if J is the first neighbor 
coupling). The fit determines h. wo=0.54±0.02 meV and 
D—IO meVA^ at 15K in both directions. 

r(q), which characterizes the excitation damping of 
the low-energy branch, obeys well a quadratic law, iden- 
tical for the [001] and [110] directions (see Fig. 9 for the 
[001] direction). 

The Q dependence of the energy-integrated intensity, 
x(Q), reveals several new features. 

1) The q-dependence of the low energy branch, mea- 
sured in the same Brillouin zone (centered at r=(002)) 
strongly differs for q//[001] and for q//[110]. This is 
illustrated in Fig. 10, where the two q directions are 
compared. The fit with a q lorentzian lineshape (con- 
tinuous lines in Fig. 10) yields correlation ranges along 
[110] and along [001] of 19 A and 4 A, respectively. This 
q-dependence is reminiscent of the ferromagnetic corre- 
lations ^// and ^j^ obtained for the static ferromagnetic 
clusters. Because of the strong decrease of the intensity 
along [110], the low-energy spin wave branch could not 
be measured beyond some q value ((^— 1.15 in Fig. 6), 
which is far from the zone boundary (C= 1.5) along this 
direction. 

(2) The intensity of the magnetic excitations of the 
low energy branch is twice larger for measurements in 
the Brillouin zone centered at t=(002), than in that cen- 
tered at t=(110). This is illustrated in Fig. 8-d, for one 
q value. Interestingly, this factor 2 corresponds to the 
geometrical factor expected for spin fluctuations trans- 



verse to the c axis, namely, the direction of ferromag- 
netism in this sample. This finding agrees with the as- 
signment of the low-energy spin-wave branch to the new 
fcrromagnetism induced by doping. We mention^that, 
in our previous study of Ca-substituted samplcst3, the 
reported measurements have suggested a different con- 
clusion (intensity larger around (110) than (002)). This 
difference is likely not real, but related to the difficulty of 
separating the domain contributions in twinned samples. 

(3) Finally, we outline the absence of any intensity 
for magnetic excitations of the high energy branch at 
-r=(002). This is illustrated in Fig. 10, where the full 
squares show the decrease of the spin wave intensities of 
the high energy branch from (001) to (002) along the [001] 
direction. This variation is a consequence of the symme- 
try of the antiferromagnetic layered structure. Therefore, 
we can assign the high energy branch to spins (or spin 
components) with the antiferromagnetic A type struc- 
ture. Accordingly, along [001], the intensity of the high 
energy branch, maximum at (001) (full squares) and the 
intensity of the low energy one, maximum at (002) (open 
circles) have complementary q-dependences. At interme- 
diate q, both types of excitations can be observed (see 
C=1.4 and 1.5 in reduced units in Fig 8-a). 

The spin dynamics of the xsr=0.06 sample has also 
been compared to that of the xca=0.05 and 0.08 samples. 

In Fig. 11, the dispersion of the low energy branch 
along [001] is reported for xca=0.05 (open squares), 
xca=0.08 (open circles) and X5'r=0.06 (full circles), all 
measurements being performed at T=15K. 
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FIG. 8. Energy scans at constant Q a): along [001] at 
k/ = 1.55A"\ b): along [001] at k/=1.3l"\ c) along [110] 
with kf — l.sA"^, d) along [110] in two different Brillouin 
zones. The continuous lines are fits with a lorentzian function 
convoluted with the spectrometer resolution. 
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FIG. 9. Linewidth of the low-energy spin wave modes F 
versus q for q along [001] at T=50K. 
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FIG. 10. q dependence of the magnon peak intensity at 
50K: along [001] for the high energy branch (filled squares) 
and the low energy branch (empty circles) and along [110] 
for the low energy branch starting from (002), and reported 
in equivalent units of rlu [001] (filled circles). The continuous 
fines are a guide for the eyes. 



FIG. 11. Dispersion of the low energy branch along [001] 
at T=15K for xca=0.05 (empty squares) xca=0.08 (empty 
circles) and xsr=0.06 (full circles). 

This comparison reveals that the low energy branch 
of xsr=0.06 corresponds very precisely to that measured 
in the xca=0.08 sample, well distinct from that in the 
xca=0.05 sample. For instance, the gap value ajo=0.55 
meV, is twice smaller than the value (0.91 meV) found 
at X(7yi=0.05. The same shift with x can be seen, coa,- 
sidering the value determined by EPR measurementsE^, 
namely luq—Q.TA meV for xsr—0-05, significantly smaller 
than that measured at the same hole concentration in the 
Ca-substituted sample (O.QlmeV). This enhancement of 
the ferromagnetism upon doping with Sr can also be seen 
when considering the antiferromagnetic coupling ("ef- 
fective" magnetic coupling J2) deduced from the high 
energy branch, which is the same as that found for 

XCa=0.08. 

We conclude that the characteristics of the low energy 
branch (gap and ferromagnetic coupling), as well as those 
of the high energy branch, evolve faster with x for Sr 
doping than for Ca doping. Moreover, the shift with x (w 
2%) is the same when comparing the spin dynamics and 
the evolution of the static ferromagnetic clusters (section 
III-B). Therefore, a relation exists between the size of 
the ferromagnetic cluster and the ferromagnetic coupling 
or the gap value defined by the low-energy spin wave 
branch. 

In this comparison, we also note that the large gap 
value of the high energy branch, 1.86 meV at T=15K, 
is, within our accuracy, the same as the gap value that 
was determined in Ca-doped samples, which is.£onstant 
in the 0.05 < xca < 0.1 concentration rangetjllj. Cor- 
respondingly, the value of the low energy branch at the 
zone boundary appears also fixed (see at (001.5) in Fig. 
11), so that the spin wave lies below the gap of the high 
energy branch, in a separated energy range. 

Preliminary results in LaMnOa+a confirm the general 
character of these observations. 




V. DISCUSSION AND CONCLUSION 

In this paper, we have shown that all the unusual static 
and dynamical features found in the Ca doped samples, 
also exist with Sr doping, confirming the general charac- 
ter of our previous observations. 

The analysis which takes into account the anisotropy of 
the spatial distribution, indicates a density of ferromag- 
netic clusters ~ 25 times smaller than the hole density. 
This suggests a picture of an electronic phase segrega- 
tion with hole-rich regions, as a general feature of the low 
doping state. The absence of any chemical segregation, 
checked by X-rays for one Ca-doped sample (x=0.08), 
confirms this electronic origin. However, the ferromag- 
netic component deduced from the Bragg-peak intensities 



cannot be explained by the volume fraction of ferromag- 
netic clusters only, suggesting that a small ferromagnetic 
component also exists within the matrix. 

The two new results are: 1) anisotropy of the ferro- 
magnetic cluster shape, indicating platelets rather than 
droplets, 2) anisotropy of the dynamic susceptibility, 
which contrasts with the isotropy of the dispersion curve. 
Both features have been observed thanks to the absence 
of twinning in our sample. 

The anisotropy of the ferromagnetic cluster shape must 
be a consequence of the AF layered magnetic, structure 
and of the double-exchange coupling of Zenercl which fa- 
vors the hole hopping within the ferromagnetic layer. 
The layered structure is related to the cooperative , 
Teller effect which lifts the degeneracy of the eg stateO 
The optimal size of the cluster can be defined from the 
competition between the kinetic energy of the e^ hole 
which is favored by the growth of the ferromagnetic clus- 
ter, and the long-range cuulomb energy, which unfavors 
the electronic segregationH. The enhanced ferromagnetic 
character induced by the Sr dopant with respect to the 
Ca dopant, confirms the role of some structural parame- 
ter for determining the hopping probability or the kinetic 
energy. 

Turning now to the spin dynamics, we observe a coex- 
istence of two spin wave branches, the high energy one, 
typical of superexchange coupling, and the low energy 
one indicating a ferromagnetic and isotropic coupling. 
From the spin- wave intensities observed at the (002) and 
(110) Bragg peaks, we can assign the high energy branch 
to spins (or spin components) with the AF A-type struc- 
ture, and the low energy one to spins (or spin compo- 
nents) ferromagnetically aligned along c. 

The fast decrease of the AF inter-layer coupling along 
[001], defined by the high-energy spin wave branch, and 
assigned to superexchange coupling, agrees well with pre- 
vious experimental resultscJ. This evolution character- 
izes an increasing ferromagnetic character, showing that 
the spin coupling , typical of the antiferromagnetic lay- 
ered structure in which the clusters are embedded, is very 
sensitive to the doping x. Specially, |-this variation with 
X, can be explained within a modelEZi, which does not 
assume any electronic phase segregation and does not 
consider the canted state. 

However, such a model cannot explain the low-energy 
spin- wave branch. The interesting result, concerning 
this low-energy spin wave branch, is the perfect isotropy 
of the corresponding dispersion (and of the damping), 
within this AF layered magnetic structure. Moreover, 
the q-dependent susceptibility, strongly anisotropic, re- 
flects closely the anisotropic shape of the ferromagnetic 
static clusters. All these characteristics lead to connect 
these magnetic excitations with the ferromagnetic hole- 
rich regions. 

The absence of any cut-off size effect or any anomaly 
at qmax in ttie dispersion curve, in-contrast with the case 
of the stripe picture for instanceE3, prevents a true in- 
homogeneous picture. Rather, it suggest that these spin 



waves are coherent through the antiferromagnetic matrix. 
The long-range ferromagnetic order, found in the system 
(canted state) is certainly responsible for that. More- 
over, the origin of the small energy gap and its evolution 
with X, could be understood by using general concepts 
of second order phase transitions. This allows to make 
some predictions for the evolution with x of x(g, w) in 
uj-q space. We recall that, as x increases, the small gap 
energy at q=0 decreases and the strength of the ferro- 
magnetic coupling defined by the q-dependence increases. 
Therefore, the gap value at x could characterize a "dis- 
tance" towards a critical Xc value at T=0. This critical 
concentration would correspond to a uniform ferromag- 
netic state. Within such an approach, the new dynamical 
susceptibility x i'-^il) can be described as a precursor ef- 
fect of this "uniform" ferromagnetic and metallic phase. 
The isotropy of the new ferromagnetic coupling indicates 
that, in this region, the magnetic unit cell is the small 
perovskite cube, as in the metallic phase occurring at 

X>Xc. 

The new ferromagnetic coupling is likely induced by 
itinerant e^ holes. Therefore, in spite of the macro- 
scopic insulatiag .character of the low-doped samples, 
band modelg3i3'B^EiE2E2l which consider itinerant Cg holes 
in addition to localized spins, should be used. 

This precursor dynamics is expected to disappear at 
the metallic transition where the hole-rich and hole-poor 
regions melt into a single phase. Actually, at X5r=0.3, a 
single spin-dynamics is observed with a vanishingly small 
gap valuetj. This will be discussed in a further paper on 
Ca substituted samples where the spin dynamics at the 
transition has been studied in a detailed wayH. 

Finally, we have outlined that, within our experimen- 
tal accuracy, the large gap of the high energy branch and 
the magnon energy of the low-energy branch at the zone 
boundary (001.5) appear to be fixed at some general val- 
ues. Moreover, the low energy branch lies below the gap 
of the high energy one (see the hatched area in Fig. 7), in 
a separated energy range. This latter peculiatieature has 
similarities with observations in La2Ni04.iil2il, where it 
is discussed in terms of host-host and impurity-host-cou- 
plings, by analogy with the single impurity problerrH. It 
is worth-noting that, unlike the case of La2NiO4.11, the 
two spin wave branches appear in two distinct symmetry 
points in the reciprocal space (001) and (002), defining 
two magnetic coupling of opposite sign along c. The com- 
petition between these two magnetic couplings could be 
responsible for this separation in energy. 

In conclusion, the static and dynamic spin correlations 
in low-doped manganites characterize a picture of two- 
magnetic regions, confirming an electronic phase segre- 
gation. However, the existence of a long range ferromag- 
netic order indicates that the ferromagnetic regions are 
connected through the matrix leading to a non uniform 
canted state. Moreover, the spin dynamics assigned to 
the AF spins or spin components, is strongly sensitive 
to the doping, acquiring a stronger ferromagnetic char- 
acter as X increases. Therefore, some homogeneous and 



10 



inhomogcncous characteristics are actually observed to- 
gether, whereas they are exclusive in all the theoretical 
models. This leads to a weakly inhomogeneous picture. 
Its true origin remains to be clarified. 
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